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Chapter 1 
Introduction 



The quark- gluon plasma (QGP) is a macroscopic system of deconfined quarks and 
gluons. The very existence of QGP at a sufficiently large temperature and/or baryon 
density is basically an unavoidable consequence of the quantum chromodynamics 
(QCD) which is a dynamical theory of strong interactions (see e.g. |^). The plasma 
has been present in the early Universe and presumably can be found in the compact 
stellar objects. Of particular interest however is the generation of QGP in relativistic 
heavy-ion collisions which has been actively studied theoretically and experimentally 
|1|] for over ten years. The life time of the plasma produced, if indeed produced, 
in these collisions is not much longer than the characteristic time scale of parton 
processe^. Therefore, QGP can achieve, in the best case, only a quasi-equilibrium 
state and studies of the nonequilibrium phenomena are crucial to discriminate the 
characteristic features of QGP. 

The transport or kinetic theory provides a natural framework to study systems 
out of thermodynamical equilibrium. Although the theory was initiated more than 
a century ago - Boltzmann derived his famous equation in 1872, the theory is still 
under vital development. Application of the Boltzmann's ideas to the systems which 
are relativistic and of quantum nature is faced with difficulties which have been 
overcome only partially till now. For a review see the monography |2^. In the case 
of the quark-gluon plasma specific difficulties appear due to the system non-Abelian 
dynamics. Nevertheless, the transport theory approach to QGP is in fast progress 
and some interesting results have been already found. 

The aim of this article is to review a few topics of the QGP transport theory. The 
first one is how to derive the transport equations of quarks and gluons. Since QCD 
is the underlying dynamical theory, these equations should be deduced from QCD. 
However, the kinetic theory of quarks and gluons have been successfully derived from 
QCD only in the mean- field or coUisionless limit |TB|, ^. The derivation of the collision 
terms is still an open question. We discuss here the issue within the dynamical model 
which is much simpler than QCD. Namely, we consider the self interacting scalar 

^The word parton is used as a common name of quarks and gluons. 
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fields with the quartic interaction term. Then, one can elucidate the essence of the 
derivation problem. 

In the third chapter we present the transport equations of quarks and gluons 
obtained in the mean-field limit. The equations are supplemented by the collision 
terms which are justified on the phenomenological ground. We briefly discuss the 
peculiarities of the transport theory of quarks and gluons and then consider the 
locally colorless plasma^. The dynamical content of QCD enters here only through 
the cross-sections of parton-parton interactions. 

The characteristic features of QGP appear when the plasma is not locally color- 
less and consequently it interacts with the chromodynamic mean field. The plasma 
response to such a field is discussed in the fourth chapter, where the color conduc- 
tivity and chromoelectric permeability tensors are found. We also analyze there the 
oscillations around the global thermodynamical equilibrium. 

The parton momentum distribution is expected to be strongly anisotropic at the 
early stage of ultrarelativistic heavy-ion collisions. Then, the parton system can be 
unstable with respect to the specific plasma modes. In the fifth chapter we discuss 
in detail the mode, which splits the parton system into the color current filaments 
parallel to the beam direction. We show why the fluctuation which initiates the 
filamentation can be very large and explain the physical mechanism responsible for the 
fluctuation growth. Then, the exponentially growing mode is found as a solution of the 
respective dispersion equation. The characteristic time of the instability development 
is estimated and finally, the possibility to observe the color filamentation in nucleus- 
nucleus collisions at RHIC and LHC is considered. 

Presenting the QGP transport theory we try to avoid model dependent concepts 
but a very crucial assumption is adopted that the plasma is perturbative i.e. the 
partons weakly interact with each other. As known, QGP becomes perturbative only 
at the temperatures much greater than the QCD scale parameter A = 200 MeV, see 
e.g. [^]. However, one believes that many results obtained in the framework of the 
perturbative QCD can be extrapolated to the nonperturbative regime. 

In the whole article we use the units where c = k = h = 1. The metric tensor is 
diagonal with goo = -gu = -g22 = -fi-ss = 1- 



^We call the plasma locally colorless if the color four-current vanishes at each space-time point. 
It differs from the terminology used in the electron-ion plasma physics, where the plasma is called 
locally neutral if the electric charge (zero component of electromagnetic four-current) is everywhere 
zero. 
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Chapter 2 



Derivation of the transport 
equation in 0^ model 



The transport equations can be usually derived by means of simple heuristic argu- 
ments similar to those which were used by Boltzmann when he formulated the kinetic 
theory of gases. However, such arguments are insufficient when one studies a system 
of complicated dynamics as the quark-gluon plasma governed by QCD. Then, one 
has to refer to a formal scheme which allows to derive the transport equation directly 
from the underlying quantum field theory. The formal scheme is also needed to spec- 
ify the limits of the kinetic approach. Indeed, the derivation shows the assumptions 
and approximations which lead to the transport theory, and hence the domain of its 
applicability can be established. 

Until now the transport equations of the QCD plasma have been successfully 
derived in the mean-field limit 3 and the structure of these equations is well un- 
derstood |18|, In particular, it has been shown that in the quasiequilibrium 
these equations provide ^ the so-called hard thermal loops The collisionless 



transport equations can be applied to the variety of problems. However, one needs the 
collision terms to discuss dissipative phenomena. In spite of some efforts p3| , 
the general form of these terms in the transport equations of the quark-gluon plasma 
remains unknown. 

The so-called Schwinger-Keldysh []65| , formulation of the quantum field theory 



provides a very promising basis to derive the transport equation beyond the mean-field 
limit. Kadanoff and Baym developed the technique for nonrelativistic quantum 



systems which has been further generalized to relativistic ones ^ |14|, |T2|, ^ 
35| , p!0| , p3| , ^ |61[] . We mention here only the papers which provide a more or less 
systematic analysis of the collision terms. 

The derivation of the complete QCD transport equations appears to be a very 
difficult task. In particular, the treatment of the massless fields such as gluons is 
troublesome. Except the well known infrared divergences which plague the pertur- 
bative expansion, there is a specific problem of nonequilibrium massless fields. The 
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inhomogeneities in the system cause the off-mass-shell propagation of particles and 
then the perturbative analysis of the collision terms appears hardly tractable. More 
specifically, it appears very difficult, if possible at all, to express the field self-energy 
as the transition matrix element squared and consequently we loose the probabilistic 
character of the kinetic theory. The problem is absent for the massive fields when the 
system is assumed to be homogeneous at the inverse mass or Compton scale. This 
is a natural assumption within the transport theory which anyway deals with the 
quantities averaged over a certain scale which can be identified with the Compton 
one. We have developed a systematic approach to the transport of massless fields. 



which allows one to treat these fields in a very similar manner as the massive ones. 
The basic idea is rather obvious. The fields which are massless in vacuum gain an 
effective mass in a medium due to the interaction. Therefore, the minimal scale at 
which the transport theory works is not an inverse bare mass, which is infinite for 
massless fields, but the inverse effective one. The staring point of the perturbative 
computation should be no longer free fields but the interacting ones. In physical 
terms, we have postulated existence of the massive quasiparticles and look for their 
transport equation. We have successfully applied the method to the massless scalar 
fields 1^3, but the generalization to QCD is far not straightforward due to the much 
reacher quasiparticle spectrum. 

To demonstrate the characteristic features of the transport theory derivation we 
discuss in this chapter the simplest nontrivial model i.e. the real massive fields with 
the lagrangian density of the form 

C{x) = - ^m'^\x) - . (2.1) 

The main steps of the derivation are the following. One defines the contour Green 
function with the time arguments on the contour in a complex time plane. This func- 
tion, which is a key element of the Schwinger-Keldysh approach, satisfies the Dyson- 
Schwinger equation. Assuming the macroscopic quasi-homogeneity of the system, one 
performs the gradient expansion and the Wigner transformation. Then, the pair of 
Dyson-Schwinger equations is converted into the transport and mass-shell equations 
both satisfied by the Wigner function. One further computes perturbatively the self 
energy which provides the Vlasov and the coUisional terms of the transport equation. 
Finally, one defines the distribution function of standard probabilistic interpretation 
and find the transport equation satisfied by this function. 

2.1 Green functions 

The contour Green function is defined as 



iA{x,y) = {T(f){x)(f){y)) 
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where (....) denotes the ensemble average at time to (usually identified with — oo); T is 
the time ordering operation along the directed contour shown in Fig. 1. The parameter 
tmax is shifted to +cxo in the calculations. The time arguments are complex with an 
infinitesimal positive or negative imaginary part, which locates them on the upper or 
on the lower branch of the contour. The ordering operation is defined as 

f(j){x)(j){y) = Q{xo,yo)(l)ix)(j){y) + Q{yo,xo)(f){y)(j){x) , 
where 0(xo, yo) equals 1 if Xq succeeds yo on the contour and equals when xq precedes 

If the field is expected to develop a finite expectation value, as it happens when 
the symmetry is spontaneously broken, the contribution {(j){x)) {(j){y)) is subtracted 



from the right-hand-side of the equation defining the Green function, see e.g. [0, ^ 
Then, one concentrates on the field fiuctuations around the expectation values. Since 
is expected to vanish in the models defined by the lagrangians ( p.!]) we neglect 
this contribution in the Green function definition. 

We also use four other Green functions with real time arguments: 

iA>{x,y) = {(j){x)(j){y)) , 

iA<{x,y) = {(j){y)(f){x)) , 
zA%x,y)''^{r<P{x)<P{y)), 
tA%x,y)''^ {T^<P{x)<P{y)) , 
where T'^{T°') prescribes (anti-)chronological time ordering: 

T''(f){x)(t){y) = e(xo - yo)(t){x)(t){y) + Qijjo - xo)(j){y)(f)c{x) , 

Ty{x)(Piy) 0(1/0 - a;o)0(x)0(y) + e(xo - Z/o)0(l/)0(x) . 
These functions are related to the contour Green functions in the following manner: 

A^{x,y) = A{x,y) for xq, yo from the upper branch, 

A"'{x,y) = A{x,y) for xq, yo from the lower branch, 

A'^{x,y) = A{x,y) for xo from the upper branch and 

yo from the lower one, 

A'^{x,y) = A{x,y) for xq from the lower branch and 

yo from the upper one. 
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It appears convenient to introduce the retarded (+) and advanced (— ) Green 
functions 

A^{x,y)''^' ±{A>ix,y) - A<ix,y))Qi±xoTyo) ■ (2.2) 

One easily finds several identities which directly follow from the definitions and relate 
the Green functions to each other. 

A'^{x, y) describes the propagation of disturbance in which a single particle is 
added to the many-particle system in space-time point y and then is removed from 
it in a space-time point x. An antiparticle disturbance is propagated backward in 
time. The meaning of A"'{x,y) is analogous but particles are propagated backward 
in time and antiparticles forward. In the zero density limit A'^{x,y) coincides with 
the Feynman propagator. 

The physical meaning of functions A^{x,y) and A^{x,y) is more transparent 
when one considers the Wigner transform defined as 

A^{X,p) = j dWP"A^(X + - -u) . (2.3) 

Then, the free-field energy-momentum tensor averaged over ensemble can be ex- 
pressed as 

C(^) = -\{^{^) <P{x)) = J -0^^p^p^tA<{X,p) . 

One recognizes the standard form of the energy- momentum tensor in the kinetic the- 
ory with the function iA^{X,p) giving the density of particles with four-momentum 
p in a space-time point X. Therefore, iA<(X,p) can be treated as a quantum ana- 
log of the classical distribution function. Indeed, the function i A< (X, p) is hermitian. 
However it is not positively definite and the probabilistic interpretation is only approx- 
imately valid. One should also observe that, in contrast to the classical distribution 
functions, iA^[X,p) can be nonzero for the off- mass-shell four- momenta. 

2.2 Equations of Motion 

The Dyson-Schwinger equations satisfied by the contour Green function are 

[a^ + m^]A{x, y) = (x, y) + /' dVn(x, x')A{x' , y) , (2.4) 

J c 

[dl + m^]A{x, y) = -6^*'^ {x, y) + j (fx'A{x, x')Il{x' , y) , (2.5) 

J c 

where n(a:, y) is the self-energy; the integration over x'q is performed on the contour 
and the function 5^'^\x^y) is defined on the contour as 

5^'^\x — y) for Xq , yo from the upper branch, 
5^^\x,y)={ for Xq , from the different branches, 

—5'-'^-' [x — y) for xq , yo from the lower branch. 
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Let us split the self-energy into three parts: 

n(x,y) = Us{x)5^'^\x,y) + U>{x,y)Q{xo,yo) + Il<{x,y)Q{yo,xo) ■ 

As we shall see later, Us provides a dominant contribution to the mean-field while 
n< determines the collision terms of the transport equations. 

With the help of the retarded and advanced Green functions (|2.2|) and the retarded 



and advanced self-energies defined in an analogous way, the equations ( p.4| ) and (|2T 
can be rewritten as 



n,(x)]A<(x,y) 



d'^x' 



n^(x, x')A~{x', y) + n+(x, x')A^(x', y) 



(2.6) 



n5(i/)]A^(x,y) 



A^(x, x')n-(x', y) + A+(a;, x')Il^{x\ y) 



(2.7) 



where all time integrations run from — oo to +00. 



2.3 Towards the Transport Equation 

The transport equations are derived under the assumption that the Green functions 
and the self-energies depend weakly on the sum of their arguments and that they 
are significantly different from zero only when the difference of their arguments is 
close to zero. For homogeneous systems, the dependence on X = {x + y) /2 drops out 
entirely due to the translational invariance and A(x,?/) depends only otiu = x — y. 
For weakly inhomogeneous, or quasihomogeneous systems, the Green functions and 
self-energies are assumed to vary slowly with X. We additionally assume that the 
Green functions and self-energies are strongly peaked near u = 0, which means that 
the correlation length is short. 

We will now convert the equations (|2.6| , p.7|) into the transport and mass-shell 
equations by implementing the above approximation and performing the Wigner 
transformation (|2.3|) for all Green functions and self-energies. This is done by means 
of the translation rules such as: 

j d^x' f{x,x')g{x\y) — 

+ 

h{x)g{x,y) — 



f{X,p)g{X,p) 



df{X,p)dg{X,p) _ df{X,p)dg{X,py 
dp^ dXi" (9X^ dp^ . 



ufv\ fv \ ^ dh{X) dg{X,p) 
hiX)g{X,p)--^^—- 
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Here d'^ = and the functions /(x, y) and g{x^ y) satisfy the assumptions discussed 
above. The function h{x) is assumed to be weakly dependent on x. 

The kinetic theory deals only with averaged system characteristics. Thus, one usu- 
ally assumes that the system is homogeneous on a scale of the Compton wave length 
of the quasiparticles. In other words, the characteristic length of inhomogeneities is 
assumed to be much larger than the inverse mass of quasiparticles. Therefore, we 
impose the condition 



(2.8) 



which leads to the quasiparticle approximation. The requirement ( p.8|) renders the 
off-shell contributions to the Green functions A< negligible. Thus, we deal with the 
quasiparticles having on-mass-shell momenta. 

Applying the translation rules and the quasiparticle approximation to Eqs. 
|2.7| ), we obtain 



p^d,--d,Mx)d^\A<{x,p) 



-(n>(x,p)A<(x,p) - n<(x,p)A>(x,p) 

^{n^(X,p),A+(X,p) + A-(X,p)} 

^{n+(x,p) + n-(x,p), A^(x,p)} , (2.9) 



+ m2-n5(X)]A^(X,p) 

^(n^(x,p)(A+(x,p) + A-(x,p)) + (n+(x,p) + n-(x,p))A^(x,p) 

(2.10) 



2 

+ ^{n>(x,p), A<(x,p)} - ^{n<(x,p), A>(x,p)} 



where we have introduced the Poisson-like bracket defined as 

(9C(X, p) dD{X, p) dC{X, p) dD{X, p) 



[C{X,p), D{X,p)] 



dp. 



dXf' 



dXi' 



dp. 



One recognizes Eq. ( |2.9| ) as a transport equation while Eq. ( p.lO| ) as a so-called 
mass-shell equation. We write down these equation in a more compact way: 

{p^ -m^ + Us{X) + Ren+(X,p), A^(X,p)} 

= z(u>{X,p)A<{X,p) - U<{X,p)A>{X,p)) 

- {n^(X,p), ReA+(X,p)}, (2.11) 



p^-m^ + Us{X) + ReU+{X,p)]A^{X,p) = -U^{X,p)ReA^{X,p) 

- Uu>{X,p), A<{X,p)} + j{ll<iX,p), A>iX,p)} . (2.12) 
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The gradient terms in the right-hand-sides of Eqs. |2.12| ) are usually neglected 
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We introduce the spectral function A defined as 

M^^y) = {[H^)^(l>iy)]) = «A>(x,y) - iA<{x,y) , 

where [0(x), denotes the field commutator. Due to the equal time commutation 
relations 

[0(t,x),0(t,y)] = , [0(t,x),0(t,y)] = -^5(3)(x- y) , 

with the dot denoting the time derivative, the Wigner transformed spectral function 
satisfies the two identities 

From the transport and mass-shell equations |2.12| ) one immediately finds 

the equations satisfied by A{X,p) which are 

-m^ + Us{X) + ReU+{X,p), A{X,p)} = 2{lmn+(X,p), ReA+(X,p)} , (2.13) 

[p"^ -m'^ + Us{X) + ReU+{X,p)]A{X,p) = 2 Imn+(X,p) ReA+(X,p) . (2.14) 
One solves the algebraic equation ( |2.14| ) as 

A(X V) = 2lmm{X,p) 

' (p2-m2 + n5(X) + Ren+(X,p))^ + (Imn+(X,p))^ ■ 

Then, it is shown that the function of the form ( |2.15|) solves Eq. (|2.13| ) as well. The 
spectral function of the free fields can be found as 

Ao(X,p) = 2ti5{p' - m2)(e(po) - e(-po)) • 

Since Ren"*" determines the quasiparticle effective mass and Imll"'" its width, the 
spectral function characterises the quasiparticle properties. 

2.4 Perturbative expansion 



As discussed in e.g. []T1|, |T^, the contour Green functions admit a perturbative 
expansion very similar to that known from the vacuum field theory with essentially 
the same Feynman rules. However, the time integrations do not run from — oo to 
+CXD, but along the contour shown in Fig. 1. The right turning point of the contour 
{tmax) must be above the largest time argument of the evaluated Green function. In 
practice, to is shifted to — oo and tmax to +oo. The second difference is the appearance 



11 



of tadpoles, i.e. loops formed by single lines, which give zero contribution in the 
vacuum case. A tadpole corresponds to a Green function with two equal space-time 
arguments. Since the Green function A{x,y) is not well defined for x = y we ascribe 
the function —iA'^{x,x) to each tadpole. The rest of Feynman rules can be taken 
from the textbook of Bjorken and Drell ||^. 

The lowest-order contribution to the self-energy, which is associated with the 
graph from Fig. 2, equals 

n(x,y) = -|5(^)(x,y)A<(x,x), 

giving 

Ils{x) = -jA<{x,x) , (2.16) 

and 

U>{x,y) = U<ix,y) = 0. 

The one-particle irreducible g"^ contributions to the self-energy are shown in Fig. 3. 
The contribution corresponding to the diagram 3a can be easily computed. However, 
it is pure real and the only effect of this contributions is a higher order modification of 
the mean-field term. Thus, we do not consider this diagrams but instead we analyse 
the contribution 3b which provides a qualitatively new effect. It gives the contour 
self-energy equal 

^c{x,y) = —Ao{x,y)Ao{y,x)Ao{x,y) , 



and consequently 

U^ix,y) = ^Al{x,y)Al{y,x)A"^{x,y) . (2.17) 



2.5 Distribution function and transport equation 

The distribution function f{X,p) is defined as 

e{po)iA<{x,p) = e(po) A{x,p) f{x,p) , 

where A{X,p) is the spectral function (|2.15|) . Then, one finds ||5^ that 

tA>iX,p) = e(po) AiX,p) if{X,p) + 1) - 0(-po) AiX,p) fix, -p) , (2.18) 



iA<{X,p) = eipo) AiX,p) fiX,p) - ei-po) A{X,p) {fix, -p) + 1) . (2.19) 



There is a very important property of expressed in the from (|2.18| , |2.19| ). 
Namely, if the Green functions A< satisfy the transport equation ( p.ll|) and the spec- 
tral function solves the equation ( |2.14| ), the mass-shell equation of A<, i.e. Eq. ( |2.12| ), 



12 



is satisfied automatically in the 0-th order of the gradient expansion. Let us note that 
the quasiparticle dispersion relation is found as a solution of the equation 



p2 - + Us{X) + ReU+{X,p) = . 



(2.20) 



The distribution function / satisfies the transport equation which can be obtained 
from Eq. ( p. 11 ) for or A^. After using Eq. (|2.13|) one finds 

A{X,p){p^-m^ + Ren+(X,p), /(X,p)} 

= iAix,p) (n>(x,p) f{x,p) - n<(x,p) (/(x,p) + 1)) 

+ z/(X,p){n>(X,p), ReA+(X,p)} 

- z(/(X,p) + 1) {n<(X,p), ReA+(X,p)} , (2.21) 

where po > 0. We have also used here the following property of the Poisson-like 
brackets: 

{A BC} = {A, B}C + {A, C}B . 

The left-hand-side of Eq. ( |2.21| ) is a straightforward generalization of the drift term 
of the standard relativistic transport equation. Computing the Poisson-like bracket 
and imposing the mass-shell constraint one finds the familiar structure 



ie(po){p'-m2 + ReU+iX,p), fiX,p)} 



= Ep(^+vV)/(X,p)-VRen+(X,p)-Vp/(X,p) , 

where the velocity v equals dEp/dp with the (positive) energy Ep being the solution 
of the dispersion equation (p.20|) . 

Let us now analyse the right- hand- side of Eq. ( |2.21| ). The collision terms are 
provided by the self energies (|2.17| ). Since the quasiparticles of interest are narrow 
(m^ -|- Ren"*" ^ Imn"*"), we take into account only those terms contributing to 
which are nonzero for the on-mass-shell momenta. The other terms are negligibly 
small 1^^. Then, the first term in r.h.s of the transport equation ( [2. 21] ) is very similar 
to the standard collision term ||2^ of the Nordheim (or Uehling-Uhlenbeck ||67|| ) 
form. Indeed, 



t{u>{x,p)f{x,p) - n<(x,p)(/(x,p) + i) 



2 



d^kAt dSA+ d^rAt 



(27r)^5(^)(p + 



k — r) 



X 



(27r)4 (27r)4 (27r)4 

ir + 1) ir + 1) f r - f r if + 1) (r + 1 



where A'^ = 0(A;o) A{X,k) and = f{X,k). The last two terms from r.h.s of 
Eq. ( p.21| ), which are neglected in the usual transport equation, are discussed in 
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Chapter 3 

Transport equations of quarks and 
gluons 



In this chapter we introduce the gauge dependent distribution functions of quarks 
and gluons. Then, we discuss the transport equations satisfied by these functions. 
Finally, a very useful notion of the locally colorless plasma is considered. 

3.1 Distribution Functions 

The (anti-)quark distribution function Q{p, x) {Q{p, x)) is a hermitian Nc x Nc matrix 
in color space (for a SU{Nc) color group) with p denoting the quark four-momentum 
and X the space-time coordinate |^ ^ The function transforms under local 
gauge transformations as 

Q{p, x) ^ U{x)Q{p, x)U\x) . (3.1) 

The color indices are here and in most cases below suppressed. 

The gluon distribution function |1^ is a hermitian (A*"^ — 1) x [N^ — 1) matrix 



^ which transforms as 

G{p,x) M{x)G{p,x)M^{x) , (3.2) 

where 

Ma,{x)=Ti[raUix)nU^ix)] 

with Ta, a = 1,...,N^ — 1 being the SU{Nc) group generators in the fundamental 
representation. One sees that, in contrast to the distribution functions known from 
the physics of atomic gases, the distribution functions of quarks and gluons have no 
simple probabilistic interpretation due to the gauge dependence. This is, however, 
not surprising if one realizes that the question about the probability to find, let us 
say, a red quark in a phase-space cell centered around {p, x) is not physical since the 
color of a quark can be changed by means of a gauge transformation. 
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It follows from the transformation laws p.2| ) that the traces of the distribution 
functions are gauge independent, and consequently they can have a familiar proba- 
bilistic interpretation. Indeed, the probability to find a quark of arbitrary color in a 
cell {p, x) is of physical meaning since it is gauge independent. The quantities, which 
are color (gauge) independent like the baryon current b'^ or the energy momentum 
tensor t^'^ , are entirely expressed through the traces of the distribution functions 



r'^ix) 



3 J (27r)3E 
d^p 



Ti[Q{p,x)] - Ti[Q{p,x)] 
Tr[g(p, x)] + Tr[Q{p, x)] + Tr[G{p, x)] 



(2TifE 

with E being the quark or gluon energy. Both quarks and gluons are assumed to be 
massless and their spin is treated as an internal degree of freedom. 

The color current, which is a gauge dependent quantity, is expressed not only 
through the traces of the distribution functions but also through the functions them- 
selves. In the Nc x matrix notation the current reads 

d^p 



nx) 



1 

-2' 



{2TxfE 



Q{p,x) - Q{p,x) - —Tt[Q{p,x)-Q{p,x) 



'^i-TafabcGbciP, x) 



(3.3) 



where g is the QCD coupling constant and fabc are the structure constants of the 



SU (Nc) group. In the adjoint representation the color current (|3.3| ) is 

^ P -''^Tr[Ta{Q{p,x) - Q{p, x))]+ifabcGbc{P, x) 



pr 



{2itYE 

where we have used the equality liijaTb 



3.2 Transport Equations 

The distribution functions of quarks and gluons satisfy the following set of transport 



equations [27, 12, 15, 16, m-. 



p^'D^Q{p,x) + gp^- 



d 1 



dp, 2 
d 1 



{F^,{x),Q{p,x)] = C[Q,Q,G] 



p'^D,Qip,x)-gp^—-{F,,{x),Q{p,x)} = C[Q,Q,G] , 

p'^V,Gip,x)+gp''^^{J',,{x),G{p,x)} = C,[Q,Q,G], (3.4) 

where {...,...} denotes the anicommutator; and are the covariant derivatives 
which act as 



T>f, = d^- ig[A^{x), ...] 
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where and A/j, are the mean-field four-potentials defined as 

Ffj_^ and J^fj^^ are the mean-field stress tensors with a color index structure analogous 
to that of the four-potentials. The mean-field is generated by the color current (|3.3| ) 
and the respective equation is 

D^F'^''{x)=fix). (3.5) 

C, C and Cg are the collision terms which vanish in the coUisionless limit, i.e. when 
the plasma evolution is dominated by the mean-field effectsQ. As already mentioned, 
the collision terms of the QGP kinetic equations have not been systematically derived 
yet and their structure remains obscure. The situation simplifies in the case of the 
colorless plasma discussed in the next section. We note that the set of transport 
equations (|3.4 |3.5|) is covariant with respect to the gauge transformations (|3.1|, |3^) . 



3.3 Colorless Plasma 

Evolving towards thermodynamical equilibrium the system of quarks of gluons tends 
to neutralize color charges. It is expected |^ that after a short period of time the 
plasma becomes locally colorless, the color current and the mean-field F^'^ vanish. 
Then, the distribution functions of quarks and gluons are proportional to the unit 
matrices in the color space. Specifically, 

Qij{p,x) = -^6ijq{p,x), i,j = l,...,Nc , 

Qij{p,x) = ■^5ijq{p,x), 

Gab{p,x) = j^^—j6abg{p,x), a,b = 1,...,N^ -1 . 

As seen the distribution functions of the colorless plasma are gauge invariant. 

The transport equations of the colorless plasma are essentially simplified. Indeed, 
taking the trace of Eqs. ( |3.4| ) one finds 

p^d^q{x,p) = c[q,q,g], 
p^d^q{x,p) = c[q,q,g], 

p^df,g{x,p) = Cg[q,q,g], (3.6) 

where c = TrC, c = TrC* and Cg = TrC^. Because the trace of a commutator is zero, 
the covariant derivatives reduce to the normal ones in (W. 



^This occurs when the characteristic mean- field frequency is much greater than the parton colh- 
sion frequency. 
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In the case of a colorless plasma the color charges can be treated as internal 
degrees of freedom of the quarks and gluons, and it is sufficient to operate with the 
color averaged quantities which are gauge independent. Then, one can imitate the 
dynamics of the colorless plasma with a non-gauge field theory model such as (j)'^. 
Then, the collision terms are of the Nordheim |^6|] (or Uehling-Uhlenbeck [^) form 
as discussed in the previous chapter. Therefore, even not knowing the collision terms 
C, C and Cg, we expect that the respective terms of the colorless plasma c, c and c^, 
which represent the binary collisions, are 



cPp2 d?Pz d^Pi 



. {2Tl fE2 (27r)3E3 (27r)3E4 

[g3g4(i - - q2) - ^ ' 



gig2(i - g3)(i - Qa. 
+ [M4(i - - - 91^2(1 - g3)(i - qa, 
+ [<?35'4(i - + 6-2) -gi5'2(i -g3)(i + 5'4 

+ [^3^4(1 - gi)(l - ^2) - gig2(l + ^3)(1 + ^4)]W^OT^<;<;(P3,P4bl, 



qq- 


-*qq 
■*qq' 


qq- 






-*q9 



P2) 



P2) 



(3.7) 



with the analogous expressions for c[q,q,g\ and Cg[g,q,q,]. We have used here the 
abbreviations qi = q{x,pi), q2 = q{x,p2) etc. Furthermore pi = p. The coefficient | 
in the first line of the r.h.s. of Eq. ( p.7| ) is required to avoid the double counting of 
identical particles. The quantities like Wqg^gg^ps, p^lpi, P2) , which corresponds to the 
quark-gluon scattering, are equal to the square of the respective matrix element mul- 
tiplied by the energy-momentum conserving 6— function. We note that the collision 
terms have to satisfy the relations 



d^p 



c[q,q,g] - c[q,q,g] 



d^p 



(27r)3E 



{2'KfE 

c[g,g, 9] + c[q,q, g] + Cg[q,q, g] 








in order to be consistent with the baryon number and energy-momentum conservation. 

In the variety of applications one uses the collision terms in the relaxation time 
approximation i.e. 



c 
c 



vp^,u^{x) (q'^'^ip, x) - q{p, x)) 
i^Pfj.u^{x) [(f'^{p, x) — q{p, X 
Ugp^u^'ix) (g^'^p, x) - g{p, x 



(3.8) 



where z/, v and Vg are the collision frequencies and is the hydrodynamic four- 
velocity which defines the local rest frame of the quark-gluon system. The equilibrium 
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distribution functions are 



2NfN, 

exp(/3^(x)p^ - (3{x)^l{x))+l ' 

2NfN, 

exp(/3'^(x)p^ + /3(x)/i(x))+l ' 

2(iV,^ - 1) 
exp(/?'^(x)p^)-l ' 

where I3'^{x) = (3{x)u'^{x), P{x) = T~^{x); T(x) and are the local temperature 
and quark chemical potential, respectively; Nf is the number of quark flavours. Spin, 
flavour and color are treated here as internal degrees of freedom. 



q^%p,x) = 
f'^{p,x) = 
g^\p,x) = 
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Chapter 4 

Plasma color response 



In this chapter we discuss how the plasma, which is colorless, homogeneous and 
stationary, responses to the color small fluctuations. 



4.1 Linear response analysis 

The distribution functions are assumed to be of the form 

Qij{p,x) = n{p)6ij + 6Qij{p,x) , (4.1) 
Qij{p, x) = n{p)5ij + 5Qij{p, x) , 
Gabip, x) = ng{p)5ab + 5Gab{p, x) , 

where the functions describing the deviation from the colorless state are assumed to 
be much smaller than the respective colorless functions. The same is assumed for the 
momentum gradients of these functions. 
Substituting ( ^ ) in ( P73D one gets 



I 7^5t^p1'^^(P'^)-'^^(P'^) (4.2) 



(27r)3E 

- ^Tr[(5Q(p, x) - 6Q{p, x)]+2iTafabJGbc{p, x) 

As seen the current occurs due to the deviation of the system from the colorless state. 
When the system becomes neutral there is no current and one expects that there is 
no mean field. Therefore, we linearize Eq. ( p.5| ) with respect to the four potential to 
the form 

d,F'^^{x)=f{x) 

with F'^'^ = d^A^ — d^A^. It should be stressed here that the linearization proce- 
dure does not cancel all non-Abelian effects. The gluon-gluon coupling, which is of 
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essentially non-Abelian character is included because the gluons contribute to the 
color current (^4.2| ). Let us also observe that in the linearized theory the color cur- 
rent is conserved (due to antisymmetry of F'^'^) i.e. d^j'^ = 0. Finally we note that, 
as shown in 0, the semiclassical QCD transport theory effectively incorporates the 
resummation over the so-called hard thermal loops |Tl[]. 

Now we substitute the distribution functions (|41|) to the transport equations ( |3.4| ) 
with the collision terms (|3.8| ). Linearizing the equations with respect to 6Q, 6Q and 
6G, one gets 



[p^d^ + UgP^u>')5G{p,x) 



-gp^'F^^ix 



gp^F^^ix 



dn{p) 

dpu 
dn(p) 

OPu 

)^^ + i^,pXK{p)-n,{p) 



+ up^u^in^'ip) - n{p)) , (4.3) 
+ vp^u^{n'"i{p) - n{p)^ , 



Performing the linearization one should remember that is of order of 6Q. Treating 
the chromodynamic field as an external one, Eqs. ( [4.3| ) are easily solved 



6Q{p,x) 
5Q{p,x) 
SG{p, x) 



-g I d'^x'Apix - x') [pf^F^^ix')^^ - upX{n''{p) " n{p))] , (4.4) 

g I d^x'Apix - x') [pf^F^^ix')^^ + upX{n'\p) ' n{p) 

,^c)ng{p) 



g / c/VAp(x-x') p^J^^.(x') 



dpu 



VgP^,u^'{nl\p) -ng{p) 



where Ap(a;) is the Green function of the kinetic operator with the collision term in 
the relaxation time approximation, 

Ap(x) = E~^Q{t) e""'' 5(=^)(x - vt) , 

with t being the zero component of x, x^ = (t,x), \ = p/E and u' = up'^u^; in the 
plasma rest frame z/' = u. 

Substituting the solutions ( [4. 41 ) in Eq. ( [4.2] ) and performing the Fourier transfor- 
mation with respect to x-variable we get 



fik)=a^^\k)Fp,ik) 
with the color conductivity tensor expressed as 



Sp 



dn{p) 



2 J i2TifE^p''{k^ + wu„) dpx 



(4.5) 



(4.6) 



+ 



dn{p) 



+ 



2NcP^pP dugip) 



p^^k^ + wua) dpx p^ik^ + WgUa) dpx 
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If the plasma colorless state is isotropic, which is the case of the global equilibrium, 
one finds that a^P^{k) = a^'^{k)u^ and Eq. ([4.5| ) gets more familiar form of the Ohm 
law, which in the plasma rest frame reads 

f{k) = a'^'^{k)E'^{k) , 

where the indices a, /?, 7 = 1, 2, 3 label the space axes and E"{k) is the a-component 
of the chromoelectric vector. The conductivity tensor describes the response of the 
QGP to the chromodynamic field. Within the approximation used here it is a color 
scalar (no color indices) or equivalently is proportional to the unit matrix in the color 
space. In the next sections we will extract the information about QGP contained in 
af'P^k). 



4.2 Chromoelectric permeability 

Let us introduce, as in the electrodynamics, the polarization vector P(x) defined as 

d 



divPfx) 



-p(x) 



dt 



P{x) =jix) 



(4.7) 



where p and j are the time-like and space-like components, respectively, of the color 
induced four-current, = (p, j). The definition ( |4.7| ) is self-consistent, only when the 
color current is conserved, not covariantly conserved. This just the case of the linear 
response approach. Further, we define the chromoelectric induction vector D(x), 

D(x) = E(x) + P(x) (4.8) 

and the chromoelectric permeability tensor, which relates the Fourier transformed D 
and E fields, 

D^ik) = e''^{k)E^{k) , (4.9) 

where = 1,2,3. Since the conductivity tensor (|]6|) is a color scalar the perme- 
ability tensor is a color scalar as well. 

Using the definitions ( [4.7| , ^l8| , ^l9|) one easily finds that 



nl3 



(k) = S' 



-cr 



aO/3 



(k) 



Fa"^^(A;) -FcT°^'^(A;) 



(4.10) 



UJ LU^ 

with a^^^{k) given by Eq. (|4.6| ); uj is the time-like component of the wave four- 
vector, k^ = (a;, k). For the isotropic plasma the two last terms in Eq. ( |4.1CI| ) vanish. 
Substituting the conductivity tensor (|4.6| ) into Eq. ( [4.10|) we get the permeability 
tensor in the plasma rest frame 



{k)=6 



+ 



3_ 

2lj 



(2^ 



dn{p) 



+ 2N, 



kv + iu dp^ 



dng{p) 



+ 



dn{p) 



— kv + iug dp^ 



CO — Viv + iu dp-y 
k'^v^i 



(4.11) 



UJ 



UJ 
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In the case of the isotropic plasma the permeabihty tensor can be expressed as 
with the longitudinal and transversal permeability functions equal 



eLik) = l + 



+ 



9 



2cjk2 J {27T 



— kv 



kvA;''' dn{p) 
uj — kv + iv dp^ 



+ 



iv dp^ 
kv k^ 



(4.12) 



dng{p) 



UJ 



kv + it'g dp^ 



d^p 



(2vr) 



1 



dn{p) 



UJ 



kv + iu dp^ 



(4.13) 



+ 



dn{p) 



kv + zz/ 



+ 



1 



dugip) 



UJ 



UJ 



+ iug dp^ 



kv F 



k2 



Because the QCD equations within the linear response approach coincide (up 
to the trivial matrix structure) with those of the electrodynamics, the dispersion 
relations of the plasma oscillations, or of plasmons, are those of the electrodynamics 
and they read g|, |6^ 

(4.14) 



det I k^^"^ - k'^k^ - uj'^e"^{k) | = . 
The relation ( [4.14| ) gets simpler form for the isotropic plasma. Namely, 

eL(A;) = 0, tT{k) = \^/uj\ 



(4.15) 



The dispersion relations determine the plasma waves which can be propagate in the 
medium. Specifically, the plane wave with co'(k), which satisfies the dispersion equa- 
tion ( [4.14| ), automatically solves the sourceless Maxwell equations in a medium. Using 
the 'quantum' language, the dispersion equation gives the relation between the energy 
and momentum of the quasiparticle excitations. In the case of plasma these are the 
transverse and longitudinal plasmons. 

There are three classes of the solutions of Eq. ( |4.14| ). Those with pure real uj are 
stable - the wave amplitude is constant in time. If the imaginary part of frequency is 
negative, the oscillations are damped - the amplitude decreases in time. Of particular 
interest are the solutions with the positive Imu; corresponding to the so-called plasma 
instabilities - the oscillations with the amplitude exponentially growing in time. 

The permeability tensor in the static limit {uj — > 0) provides the information about 
the plasma response to constant fields. Computing eiiuj = 0,k) for the equilibrium 
collisionless plasma one finds 



ei(u; = 0,k) = 1 + 



m 



D 



k2 
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where tjid is the Debye mass which for the baryonless plasma of massless quarks and 
gluons equals 



^2 _ g^T\Nf + 2N, 



< = - '^'^^ "^-^ . (4.16) 

The chromoelectric potential of the static point-like source embeded in the plasma, 
which equals Q 



^o(x) = — rexp{-mD I X I) , 
47r I X I 

is screened at the distance l/mo- 

Since the parton density is ~ T^, one finds from Eq. ( [4.16| ) that the number of 



partons in the Debye sphere (the sphere of the radius equal to the screening length) 
is ~ It is much greater than unity if the plasma is perturbative i.e. when 

1/g ^ 1. A large parton number in the Debye sphere justifies, in particular, the 
use of the mean-field to describe QGP. Let us also mention that the ultrarelativistic 
perturbative plasma is automatically ideal i.e. the average parton interaction energy, 
which is ~ g"^/ (r) with (r) ~ being the average interparticle distance, is much 
smaller than the parton thermal energy which equals ~ T. This is not the case for 



the nonrelativistic electron plasma. Then, the screening length is (see e.g. p6i |64[ ) 

mj^ = e — , 

with He being the electron density^, which is independent of the temperature. As seen, 
the smallness of the coupling constant does not guarantee that the nonrelativistic 
plasma is ideal. This occurs when the number of electrons in the Debye sphere is 
large, i.e. when T^^^ ^ e^rig/^. 



4.3 Oscillations around the global equilibrium 

Substituting the equilibrium distribution functions, Fermi-Dirac for quarks and Bose- 
Einstein for gluons, into Eqs. (|4.13| ) and (|4.14| ) one finds the permeability functions 
ei and e^, which for the coUisionless {u = u = Ug = 0) and baryonless (/i = 0) plasma 
of massless partons can be computed analitically as 



eL = l + 



fc2 



1-4 
2k 



e-r = 1 



^2k 



In 



2u^ 



k + u ^ 
k — uj 

, , k + uj 
m 



iTiQ{k 



k — UJ 



-iTiQ{k 



00] 



(4.17) 



(4.18) 



^We use the units, where the fine structure constant a — e^/47r. In the Gauss units, which are 
traditionally used the electron-ion plasma physics, a = . 
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where k 



k I and uq is the plasma frequency equal 

g^T\Nf + 2N,) 
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(4.19) 



One sees that for lu > k the dielectric functions (|4.171 , |4.18|) are purely real i.e. there 
are no dissipative processes. 

Substituting ([4.17| , |4.18|) into ([4.15|) one finds the dispersion relation of the longi- 
tudinal mode (the chromoelectric field parallel to the wave vector) 



UJ 



^0 + i^^ 

A;2(l + 4exp(-2 



2ky3ujl 



for c<Jo ^ /c 
for ujq k 



and 



ooi + Ik' 



for Uq k 
^^ ^uJq + k'^ , for ooo k 

for the transverse one (the chromoelectric field perpendicular to the wave vector). 



Because the longitudinal and transverse oscillations are time-like 



> k^), the 



phase velocity of the waves is greater than the velocity of light. For this reason the 
Landau damping is absent. As known, the Landau damping is due to the coUisionless 
energy transfer from the wave to the plasma particles, the velocity of which is equal 
to the wave phase velocity |^ . 



The oscillations of the coUisionless QGP around global equilibrium have been 
studied by means of the transport theory in several papers |l^, |, |4^. The 



problem has been also discussed in |30, 31] using a specific variant of the QGP theory 



with the classical color |27, 



In the above presentation we have followed [49 



The dispersion relations given above agree with those found in the finite-temperature 

mm. 



QCD within the one-loop approximation, see e.g. ||7y, |_ 

Let us now consider the dielectric function with nonzero equilibration rates. As 
previously, the partons are massless and the plasma is baryonless which imposes 
= u. Then, one easily evaluates the integrals (|4.13| ) and ( [4.14|) for uj ^ k, u ^ u 



and CO Ug. The results read 



-02 + 

V 5 



for the longitudinal mode and 



2 2 /•2 I J 2 I ^ T 2 

^ = ^0 -C +70 +^k 
4 5 



for the transverse one; u and 7 denote the real and imaginary part, respectively, of 
the complex frequency; and ( are parameters related to the equilibration rates. 



Nf + 2N, 



+ 1^, 



2Nr 



^Nf + 2Nc 



(4.20) 



A^ 



/ 



Nf + 2N, 



2N, 
Nf + 2Nc 
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One sees that, when compared with the colhsionless plasma, the frequency of the 
oscillations is smaller and that the oscillations are damped. To find the numerical 
value of the damping rate - the plasma oscillation decrement 7, the equilibration rates 
(z/ and Ug) have to be estimated. 

If u or Ug is identified with the mean free flight time controlled by the binary 
collisions, the equilibration rate is of the order g'^lnl/g. However, in the relativistic 
plasma there is another damping mechanism which is the plasmon decay into quark- 
antiquark or gluon-gluon pair. It is easy to observe that, even in the limit of massless 
quarks, the decay into gluons is much more probable than that into quarks [^, ^ |. 



Let us consider the decay of plasmon of zero momentum. The phase-space volume of 
the decay final state is proportional to 

{lTn''{uJo/2))\ (4.21) 

where the upper sign refers to quarks while the lower one to gluons. Since the plasma 
frequency ([4.19|) is much smaller than the temperature in the perturbative plasma. 



the factor ( ^.21| ) can be expanded as 

, 2 



(l-n''^{uJo/2)) ^ l/4 + ^o/8T 
(l + n^^K/2))' = AT'/u^. 



One sees that the decay into gluons is more probable than that into quarks by a factor 
of order g^"^ . 



Using the standard rules of finite-temperature field theory, one easily finds (see 
e.g. ii) the width of the zero-momentum plasmon due to the decay into gluons 



2^371 



which is the same for longitudinal and transverse plasmons. However, cannot 
be identified with the plasmon equilibration rate F, because the plasmon decays are 
partially compensated by the plasmon formation processes. As shown in |7ll], see also 
the formation rate Fj is related to F^^ as 



Tj = exp(-a;o/T)F, ^ (1 - a;o/T)F, . 
Since the equilibration rate of the plasmon F = F^^ — Fj [|7T|, one finds 



(4.22) 



We note that F^ and Tj are of the order of g, while F is of the order of g^. Since 
there is a preferred reference frame - the rest frame of the thermostat, the plasmon 
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decay width is not a Lorentz scalar. Therefore, the resuh ( [4.22| ) is vahd only for the 
zero-momentum or approximately long-wave plasmons. 



Substituting Ug equal ( |4.22| ) and z/ = into Eq. ( |4.2CI| ), one estimates the decre- 



ment of the oscillation damping as 

Although u = the damping rate depends on the number of quark flavours. This 
seems to be in agreement with the physical intuition. When the number of quark 
flavours is increased the inertia of the system is also increased, and consequently the 
time needed to damp the oscillations is longer. However Eq. ([4.23|) disagrees (by a 
factor 2Nc/ {Nf+2Nc)) with the result from [^, where 7 equals ( 4.22| ). Unfortunately, 



the discrepancy can not be resolved within the relaxation time approximation and a 
more elaborated analysis is needed. 
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Chapter 5 

Filamentation instability 



In the near future the nucleus-nucleus collisions will be studied experimentally at 
the accelerators of a new generation: Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven and Large Hadron Collider (LHC) at CERN. The collision energy will 
be larger by one or two orders of magnitude than that one of the currently operating 
machines. A copious production of partons, mainly gluons, due to hard and semihard 



processes is expected in the heavy- ion collisions at this new energy domain [22, 58 



Thus, one deals with the many-parton system at the early stage of the collision. The 
system is on average locally colorless but random fluctuations can break the neu- 
trality. Since the system is initially far from equilibrium, specific color fluctuations 
can exponentially grow in time and then noticeably influence the system evolution. 
While the very existence of such instabilities, similar to those which are known from 
the electron- ion plasma, see e.g. H], is fairly obvious and was commented upon long 



time ago |28|, it is far less trivial to find those instabilities which are relevant for the 



parton system produced in ultrarelativistic heavy-ion collisions. 



A system of two interpenetrating beams of nucleons |3^ or partons |5^, ^ 



58| , Q was argued to be unstable with respect to the so-called filamentation 
or Weibel instability [^. However, such a system appears to be rather unrealistic 
from the experimental point of view. Then, we have argued [^, ^ ^ that the 
filamentation can occur under weaker conditions which are very probable for heavy- 
ion collisions at RHIC and LHC. Instead of the two streams of partons, it appears 
sufficient to assume a strongly anisotropic momentum distribution. We systematically 
review here the whole problem. 



5.1 Fluctuation spectrum 

We start with the discussion on how the unstable modes are initiated. Specifically, 
we show that the fiuctuations, which act as seeds of the filamentation, are large, much 
larger than in the equilibrium plasma. Since the system of interest is far from the 
equilibrium, the fluctuations are not determined by the chromo dielectric permeability 
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tensor discussed in the previous section. The fluctuation-dissipation theorem does 
not hold in such a case. Thus, we derive the color current correlation function which 
provides the fluctuation spectrum. 

QGP is assumed to be on average locally colorless, homogeneous and stationary. 
Therefore, the distribution functions averaged over ensemble are of the form 

{Q^j{t, X, p)) = 6'^n{p) , {Qij{t, X, p)) = 6'^n{p) , 

(G,,(t,x,p)) = 5'^^n,(p), 

which give the zero average color current. 

We study the fluctuations of the color current generalizing a well-known approach 
to the fluctuating electric current @. For a system of noninteracting quarks and 
gluons we have derived (in the classical limit) the following expression of the current 
correlation tensor 

ACa.x) = (i;(i,,x,bt((2,x.)) 

where the effective parton distribution function /(p) equals n(p) + n{p) + 2Njig{p) 
and (t,x) = (^2 — tii^2 — Xi). Due to the average space-time homogeneity the 
correlation tensor depends only on the difference (^2 — ^11^2 — xi). The physical 
meaning of the formula ( |5.1| ) is transparent. The space-time points (ti, xi) and (i:2, ^2) 
are correlated in the system of noninteracting particles if the particles fly from (ti, xi) 
to (^25X2). For this reason the delta (5(^''(x — vt) is present in the formula ( ^.1|) . The 
momentum integral of the distribution function simply represents the summation over 



particles. The fluctuation spectrum is found as a Fourier transform of the tensor ([O 



I.e. 



Mr(^, ^) = \9"^''' j ^/-^ /(P) 2vr5(^ - kv) . (5.2) 

When the system is in equilibrium the fluctuations are given, according to the 
fluctuation-dissipation theorem, by the respective response function. For /(p) be- 
ing the classical equilibrium distribution function one indeed finds the standard 
fluctuation-dissipation relation valid in the (7^— order. For example, 

Mr(a.,k) = <5'^^^-Ime,.(a;,k), 

277 UO 

where T is the temperature and represents the longitudinal part of the chromodi- 



electric tensor ([4.13| ). 
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5.2 Parton distributions 



We model the parton momentum distribution at the early stage of ultrarelativistic 
heavy-ion collision by two functions: 

/(P) = i^QiY - y)e{Y + y) h{p^) -\- , (5.3) 
2Y p± chy 

and 

/(P) = ^Qi-P - p\\)Q{V + p\\) h{p^) , (5.4) 

where p\\ and p_L denote the parton rapidity, the longitudinal and transverse mo- 
menta, respectively. The parton momentum distribution ( |5.3| ) corresponds to the 
rapidity distribution which is flat in the interval (— Y,y). The distribution ( |5.4| ) is 
flat for the longitudinal momentum —V < p\\ < V. We do not specify the trans- 
verse momentum distribution h{p±), which is assumed to be of the same shape for 
quarks and gluons, because it is sufficient for our considerations to demand that the 
distributions (|5.3| , |5^) , are strongly elongated along the z— axis i.e. ^ 1 and 
V > {pi_). 

The QCD-based computations, see e.g. ||22|, |6^, show that the rapidity distri- 
bution of partons produced at the early stage of heavy-ion collisions is essentially 
gaussian with the width of about one to two units. When the distribution ( |5.3|) sim- 
ulates the gaussian one, Y does not measure the size of the 'plateau' but rather the 
range over which the partons are spread. If one takes the gaussian distribution of the 
variance a and the distribution (|5.3|) of the same variance, then Y = y/3 a. 

As already mentioned, the parton system described by the distribution functions 
( p73| , |5^) is assumed to be homogeneous and stationary. Applicability of this assump- 
tion is very limited because there is a correlation between the parton longitudinal 
momentum and its position, i.e. partons with very different momenta will find them- 
selves in different regions of space shortly after the collision. However, one should 
remember that we consider the parton system at a very early stage of the collision, 
soon after the Lorentz contracted ultrarelativistic nuclei traverse each other. At this 
stage partons are most copiously produced but do not have enough time to escape 
from each other. Thus, the assumption of homogeneity holds for the space-time do- 
main of the longitudinal size, say, 2 — 3 fm and life time 2 — 3 fm/c. As shown 
below, this time is long enough for the instability to occur. 



5.3 Seeds of the filament at ion 

Let us now calculate the correlation tensor (15. 21) for the distribution functions 



5^ . Due to the symmetry /(p) = /(— p) of these distributions, the tensor M^'^ is 
diagonal i.e. M'^" = for ii ^ v. Since the average parton longitudinal momentum 
is much bigger than the transverse one, it obviously follows from Eq. (|5.2|) that the 
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largest fluctuating current appears along the axis. Therefore, we discuss the M^^ 
component of the correlation tensor. M^^{uj, k) depends on the k— vector orientation 
and there are two generic cases: k = (A;^,0,0) and k = (0,0,^2). The inspection of 
Eq. (|5.2| ) shows that the fluctuations with k = {kx, 0, 0) are much larger than those 
with k = (0, 0, fcz). Thus, let us consider M^^{uj,kx)- Substituting the distributions 
( p73| , p.4| ) into ( |5.2| ) one finds after azimuthal integration that M^^{u, k^) reaches the 
maximal values for uj"^ <^ k"^. So, we compute M^^ at a; = 0. Keeping in mind that 
^ 1 and V ^ {p±) we get the following approximate expressions for the flat y— 
and distributions: 

1 In) 

AG(u; = 0,AJ = -j^*»'-^, (5^5) 



M^(. = 0.fe) = i/r'^M, (5.6) 



where (p) is the effective parton density given for iVc = 3 as 

= / W? ^^^^ " 4^ i '^^^^^ ^^^^^ " \ ^ I ' 



with {p)gg denoting the average density of quarks and antiquarks, and {p)g that of 
gluons. For the flat pn— case we have also used the approximate equality 



poo ^ roo 

/ dp± h{p±) = - — - / dp±ps_ h{p, 
Jo {V\) Jo 



{P±) -Jo 

to get the expression ( ^.6|) . It is instructive to compare the results { f).5[ |5^ ) with the 
analogous one for the equilibrium plasma which is 

M-{u; = 0,kx) = ^g' S^"!^^. 

One sees that the current fluctuations in the anisotropic plasma are amplified by the 
large factor which is /Y or V/ {p±)- With the estimated value of Y 2.5 for RHIC 
and 5.0 for LHC 0, the amplification factor /Y equals 4.9 and 29.7, respectively. 

5.4 Filamentation mechanism 

Following [|l^ we are going to argue that the fluctuation, which contributes to 
Ml^iuj = 0, kx), grows in time. The form of the fluctuating current is 

jaix) = ja cos(A;^a;) , (5.7) 
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where is the unit vector in the direction. Thus, there are current filaments 
of the thickness TT/\kx\ with the current flowing in the opposite directions in the 
neighboring filaments. For the purpose of a qualitative argumentation presented here 
the chromodynamics is treated as an eightfold electrodynamics. Then, the magnetic 
field generated by the current (|5.7| ) is given as 

Ba(a;) = sin(fc^a;) , 

while the Lorentz force acting on the partons, which fly along the beam, equals 

F(a;) = ga V X Ba{x) = -qa ^ sin(A;^x) , 

i^x 

where is the color charge. One observes, see Fig. 4, that the force distributes the 
partons in such a way that those, which positively contribute to the current in a given 
filament, are focused to the filament center while those, which negatively contribute, 
are moved to the neighboring one. Thus, the initial current grows. 



5.5 Dispersion equation 

We analyse here the dispersion equation which for the anisotropic plasma is of the 
form ([4.14|) . The plasma is assumed to be collisionless i.e. the mean-field interaction 
dominates the system dynamics and v = iO~^. The assumption is correct if the inverse 
characteristic time of the mean-field phenomena is substantially larger than the 
collision frequency u. Otherwise, the infinitesimally small imaginary quantity iO"*" 
from (|4.12|) should be substituted by iu. Such a substitution however seriously com- 



plicates analysis of the dispersion equation ( |4.14| ). Therefore, we solve the problem 
within the collisionless limit and only a posteriori argue validity of this approxima- 
tion. 



As already mentioned, the solutions t^(k) of ( [4.14|) are stable when Imu < and 



unstable when Imui > 0. It appears difficult to find the solutions of Eq. ( [4.14 ) 



because of the complicated structure of the chromodielectric tensor ( [4.12|) . However, 
the problem simplifies because we are interested in the specific modes with the wave 
vector k perpendicular and to the chromoelectric field E parallel to the beam. Thus, 
we consider the configuration 

E=(0,0,E,), k=(fc,,0,0). (5.8) 

Then, the dispersion equation (|4.14| ) gets the form 



H{uj) = kl - uh'^icu, k^)=0 , (5.9) 
where only one diagonal component of the dielectric tensor enters. 
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5.6 Penrose criterion 



The stability analysis can be performed without solving Eq. ( ^.9|) explicitly. Indeed, 
the so-called Penrose criterion states that the dispersion equation H{uj) = has 
unstable solutions if H{uj = 0) < 0. The meaning of this statement will be clearer 
after we will approximately solve the dispersion equation in the next section. 
Let us compute H{0) which can be written as 

H{0) = kl-x', (5.10) 

with 



where the plasma frequency parameter is 



As we shall see below, uq gives the frequency of the stable mode of the configuration 



(PD when k^^O. 

Substituting the distribution functions (|5.3| , |5l^) into Eqs. (|5.11|) and ( ^.12|) one 



finds the analytical but rather complicated expression of H{0). In the case of the flat 
y— distribution we thus take the limit ^ 1, while for the flat py— distribution we 
assume that V ^ (pi.)- Then, we get for the fiat y— distribution 

^ £l f,Jhip,^+,/J^] = f^i £; prMtf") , (5.13) 



An Y J dp± J An Y 

and for the flat pn— distribution 



where as = g'^/An'^ is the strong coupling constant and p™™ denotes the minimal 
transverse momentum. The function h{pj_) is assumed to decrease faster than l/p± 
when p± ^ oo. 

As seen, the sign of H{0) given by Eq. ( S.lOj) is (for sufficiently small k"^) deter- 



mined by the transverse momentum distribution at the minimal momentum. There 
are unstable modes if p™™/;,(p™™) > for the flat y— distribution and if h{pf_™) > 
for the flat p\\ case. Since the distribution h{p±) is expected to be a monotonously 
decreasing function of p±, the instability condition for the flat py— distribution seems 
to be always satisfled. The situation with the flat distribution is less clear. So, let 
us discuss it in more detail. We consider three characteristic cases of h{p±) discussed 
in the literature. 
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1. The transverse momentum distribution due to a single binary parton-parton 
interaction is proportional to pj^ |19| and blows up when p± 0. In such a 



case > 0, there are unstable modes and p™™ should be treated as a 

cut-off parameter reffecting e.g. the finite size of the system. 

2. The transverse momentum distribution proportional to {p± + m^)"^'^ with 
m_L = 2.9 GeV has been found in 0, where except the binary parton-parton 
scattering the initial and final state radiation has been taken into account. This 
distribution, in contrast to that from 1), gives = for p™™ = and 
there is no instability. Although one should remember that the finite value of 
■m± found in ^ is the result of infrared cut-off. Thus, it seems more reasonable 
to use the distribution from 1), where the cut-off explicitly appears. 

3. One treats perturbatively only partons with p±_ > assuming that those 
with lower momenta form colorless clusters or strings due to a nonperturbative 
interaction. It should be stressed that the colorless objects do not contribute 
to the dielectric tensor (|4.12|) , which is found in the linear response approxi- 
mation [^, Thus, only the partons with p±_ > p™" are of interest for us. 
Consequently is positive and there are unstable modes. As shown 
in [Q, the screening lengths due to the large parton density are smaller than 
the confinement scale in the vacuum. Therefore, the cut-off parameter p™"^ 
should be presumably reduced from 1-2 GeV usually used for proton-proton 
interactions to, let us say, 0.1 - 0.2 GeV. 

We cannot draw a firm conclusion but we see that the instability condition is triv- 
ially satisfied for the fiat p\\ —distribution and is also fulfilled for the fiat ?/— distribution 
under plausible assumptions. Let us mention that the difference between the instabil- 
ity conditions for the fiat y— and py— distribution is due to a very specific property of 
the y— distribution which is limited to the interval {—Y, Y). The point is that y — >■ ibcxo 
when p± —>■ and consequently, the limits in the rapidity suppress the contribution 
from the small transverse momenta to the dielectric tensor. For this reason we need 
for the instability the distribution h{p±) which diverges for p± ^ in the case of 
the fiat ?/— distribution, while the instability condition for the fiat py— distribution is 
satisfied when h{0) is finite. If we assumed the gaussian rapidity distribution instead 
of (|5.3D, the instability condition would be less stringent. In any case, we assume 
that the Penrose criterion is satisfied and look for the unstable modes solving the 
dispersion equation (|5.9| ). 
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5.7 Unstable mode 

The dispersion equation ( ^.91) for a cylindrically symmetric system is 

dpii pf. 



a. 



Att'^ Jo 

df r2. 



dp± 



p\ + pi 



COS(5 







dp±_ Jo a — COS0 + iO+ 
with the plasma frequency Uq given by Eq. ( |5.12| ) and a denoting 



(5.15) 



u \jp\\+pi 

kx P± 

We solve Eq. ( p.l5[ ) in the two limiting cases \uj/kx\ ^ 1 and \kx/uj\ ^ 1. In the 
first case the azimuthal integral is approximated as 



2tt 



cose 



IT 



a — COS0 + iO+ 



+ Oia- 



Then, the equation (|5.15| ) gets the form 



; 2 2 I 2 I n 



(5.16) 



where r/, as cuq, is a constant defined as 



„2_ ^ , , pIpI ^/(p 



We have computed ujq and r/ for the fiat pn— and distribution. In the limit 



1 and P ^ respectively, we have found 

uJa = z^-:: I dp±p±h{pA 



2'kV 



and 



a. 



16Y 



In 



47rP V(^'^), 



dp±p_i 



2 c?^(p±) 



(5.17) 
(5.18) 

(5.19) 
(5.20) 
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The solutions of Eq. ( |5.16| ) are 



\ {e + u^l ± ^{kl + u^lf + Ar^^kl ) . (5.21) 



One sees that u^^ > and a;^ < 0. Thus, there is a pure real mode a;+, which is 
stable, and two pure imaginary modes cu-, one of them being unstable. As mentioned 
previously, uoj^ = uq when kx = 0- 

Let us focus our attention on the unstable mode which can be approximated as 

2 ^ f -^kl , for kl < Lui 
[ -7]^ , for kl > ^2 



One should keep in mind that Eq. ( [5.21| ) holds only for \uj/kx\ ^ 1. We see that 



can satisfy this condition for k^. <C lOq if rj'^ 3> Wq and for k^ ^ cug if r/^ -C Wq. To 
check whether these conditions can be satisfied, we compare rj"^ to u^. Assuming that 
h{p±) ~ pj_^, one finds from Eqs. ( ^.19| , ^.20| ) 



V = — 3 — ^0 • (5-22) 



2 



^0 



Since /3 = 6 |TP|, |^ we get 77^ > 3uJq. Therefore, the solution ( p.21| ) for k^ <^ Uq should 
be correct. 

Let us now solve the dispersion equation (|5.15| ) in the second case when \kx/uj\ ^ 
1. Then, the azimuthal integral from Eq. ( p.l5| ) is approximated as 

d(j) cosd) , , 

Jo a — COS0 + «L)+ 

and we immediately get the dispersion relation 

u;' = kl-x\ (5.23) 



with given by Eq. ( |5.13| ) or ( |5.14| ). As previously we have assumed that ^ 1 



and V ^ {p±)- Eq. ( ^.23|) provides a real mode for k^. > and two imaginary modes 
for kl < x^. Since the solution ( p.23| ) must satisfy the condition \kx/uj\ 1, it holds 
only for kl > \kl - x^\- 

The dispersion relation of the unstable mode in the whole domain of wave vectors 
is schematically shown in Fig. 5, where the solutions ( p.21| ) and ( |5.23| ) are combined. 



Now one sees how the Penrose criterion works. When = the unstable mode 
disappears. 
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5.8 Time scales 



The instability studied here can occur in heavy-ion colhsions if the time of instabihty 
development is short enough, shorter than the characteristic time of evolution of the 
nonequilibrium state described by the distribution functions ^.41) . 

Let us first estimate the time of instability development which is given by 1/Imci;. 
As seen in Fig. 5, |Imix;| < rj. Thus, we define the minimal time as Tmin = ^/v- To find 
Tmin we estimate the plasma frequency. We consider here only the flat |/— distribution 
which seems to be more reasonable than the fiat py— distribution. Approximating 
/ dp±h{p±) as / dp±p±h{pj_) / {p±) the plasma frequency ( ^.12|) can be written as 



^o = QY4^iN, + N, + -N,), (5.24) 

where Nc = 3; Nq, Nq and A^^^ are the numbers of quarks, antiquarks and gluons, 
respectively, produced in the volume, which has been estimated in the following way. 
Since we are interested in the central collisions, the volume corresponds to a cylinder 
of the radius tqA^^^ with tq = 1.1 fm and A being the mass number of the colliding 
nuclei. Using the uncertainty principle argument, the length of the cylinder has been 
taken as 1/ {p±), which is the formation time of parton with the transverse momentum 

{P±)- 

Neglecting quarks and antiquarks in Eq. ( ^.24|) and substituting there Ng = 570 
for the central Au-Au collision at RHIC {Y = 2.5) and Ng = 8100 for the same 
colliding system at LHC {Y = 5.0) 0, we get 

LUo = 280 MeV for RHIC , luq = 430 MeV for LHC 

for as = 0.3 at RHIC and a, = 0.1 at LHC. Using Eq. ( ^1221) with p = 6 one finds 



Tinin = 0.4 fm/c for RHIC , t^i^ = 0.3 fm/c for LHC . 

The plasma has been assumed collisionless in our analysis. Such an assumption 
is usually correct for weakly interacting systems because the damping rates of the 
collective modes due to collisions are of the higher order in as than the frequencies 
of these modes, see e.g. 0. However, it has been argued recently [^] that the color 



collective modes are overdamped due to the unscreened chromomagnetic interaction. 
However it is unclear whether these arguments concern the unstable mode discussed 



here. The point is that the paper deals with the neutralization of color charges 
which generate the longitudinal chromoelectric field while the unstable mode which 
we have found is transversal and consequently is generated by the color currents not 
charges. Let us refer here to the electron-ion plasma, where the charge neutralization 
is a very fast process while currents can exist in the system for a much longer time 



43| . In any case, the above estimates of the instability development should be treated 



as lower limits. 
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Let us now discuss the characteristic time of evolution of the nonequihbrium 
state described by the distribution functions (|5.3| , pl4|) . Except the possible unstable 
collective modes, there are two other important processes responsible for the temporal 
evolution of the initially produced many-parton system: free streaming ^, ^ 
and parton-parton scattering. The two processes lead to the isotropic momentum 
distribution of partons in a given cell. The estimated time to achieve local isotropy 
due to the free streaming is about 0.7 fm/c at RHIC ||20|. The estimates of the 



equilibration time due to the parton scatterings are similar [21, 36 1. As seen the three 



time scales of interest are close to each other. Therefore, the color unstable modes 
can play a role in the dynamics of many-parton system produced at the early stage 
of heavy-ion collision, but presumably the pattern of instability cannot fully develop. 



5.9 Detecting the filament at ion 

One asks whether the color instabilities are detectable in ultrarelativistic heavy-ion 
collisions. The answer seems to be positive because the occurrence of the filamentation 
breaks the azimuthal symmetry of the system and hopefully will be visible in the final 
state. The azimuthal orientation of the wave vector will change from one collision 
to another while the instability growth will lead to the energy transport along this 
vector (the Poynting vector points in this direction). Consequently, one expects 
significant variation of the transverse energy as a function of the azimuthal angle. 
This expectation is qualitatively different than that based on the parton cascade 



simulations ||22|, where the fluctuations are strongly damped due to the large number 
of uncorrelated partons. Due to the collective character of the filamentation instability 
the azimuthal symmetry will be presumably broken by a flow of large number of 
particles with relatively small transverse momenta. The jets produced in hard parton- 
parton interactions also break the azimuthal symmetry. However, the symmetry is 
broken in this case due to a few particles with large transverse momentum. The 
problem obviously needs further studies but the event-by-event analysis of the nuclear 
collision seems to give a chance to observe the color instabilities in the experiments 
planed at RHIC and LHC. 
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Figure Captions 

Fig. 1. The contour along the time axis for an evaluation of the operator expectation 
values. 

Fig. 2. The lowest-order diagram of the self-energy. 
Fig. 3. The second-order diagrams of the self-energy. 

Fig. 4. The mechanism of filamentation. The phenomenon is, for simplicity, consid- 
ered in terms of the electrodynamics. The fluctuating current generates the magnetic 
field acting on the positively charged particles which in turn contribute to the current 
(see text). <S> and denote the parallel and, respectively, antiparallel orientation of 
the magnetic field with respect to the y— axis. 

Fig. 5. The schematic view of the dispersion relation of the filamentation mode. 
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